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Abstract 



/\ • The transition from quasiperiodicity to chaos is studied in a two-dimensional 



dissipative map with the inverse golden mean rotation number. On the basis of 
a decimation scheme, it is argued that the (minimal) slope of the critical iterated 
circle map is proportional to the effective Jacobian determinant. Approaching 
the zero-Jacobian-determinant limit, the factor of proportion becomes a universal 
constant. Numerical investigation on the dissipative standard map suggests that 
this universal number could become observable in experiments. 



In an almost linear continuous dynamical system with two competing frequences, 
the attractor is typically either a mode-locked periodic state or an invariant 
torus associated with quasiperiodic motion. In the parameter space, periodic and 
quasiperiodic attractors are mingled up in such a way that there is a positive 
probability for each kind of orbit to be found. Experiments indicate that mode- 
locking and quasiperiodic behaviour are generic in hydrodynamics [1], charge- 
density-wave conductors [2] and other physical [3,4] and chemical [5] systems. 
Changing parameter values may increase the nonlinearity and lead to a transition 
to chaos. 

The experiment by Martin and Martienssen [3] on the electrical conductivity of 
barium sodium niobate crystals is a very nice example of the case in which it is 
possible to measure the actual return map characterizing the discretized dynamics 
on an invariant circle (the invariant torus appears as an invariant circle for the 
Poincare map of the system). Bohr et al. [6,7] point out the intimate connection 
between the existence of an invariant circle and the one-dimensional nature of 
the return map. In particular, they show that a zero slope in the return map is 
impossible if the underlying invariant circle is smooth. The fact that the invariant 
circle loses smoothness before breaking up [8,9] could mislead one into thinking 
that a zero slope in the return map is a necessary condition for the system to 
be critical, i.e. about to become chaotic. Another motivation for this kind of 
false idea could come from the fact that an analytic circle map has a zero-slope 
inflection point at the transition from quasiperiodicity to chaos [10]. However, if 
the Jacobian determinant of the Poincare map is positive everywhere along the 
invariant circle, it is impossible that the "reduced" circle map, i.e. the projection 
of the Poincare map on the invariant circle, would have a zero slope at some point 
of the circle. The reduced circle map could have a zero slope only if the tangent 



vector at the corresponding point was annihilated by the Jacobian matrix. This 
could happen only if the Jacobian determinant vanished. The positivity of the 
slope of a critical circle map has been noticed by several authors [11]. 

In this letter the relation between the slope of the reduced circle map and the 
Jacobian determinant is elaborated further. For simplicity, I will restrict myself 
to a two-dimensional Poincare map with rotation number C = (y^ — l)/2. The 
reduced circle map is denoted by h{x) where a: is a scaled "angle" variable for the 
invariant circle so that h{x + 1) = h{x) + 1. By the assumption on the rotation 
number, h'^{x)/n -^ C and, moreover, hn{x) — Fn-i =h^^ (x) — F^-i^' x as n tends 
to infinity. Here F^ stands for the nth Fibonacci number, Fn+i = F^ + -Fn-i 
{Fq = 0, Fi = 1). The Jacobian determinant of the F^ times iterated map 
at X is denoted by Jn{x). It will be shown below that in the critical case 
h'^{xo) ~ Jn{xo), where xq is a special point [8,9] associated with the universal 
scaling by a ~ —1.2885746. xq corresponds to a cubic critical point for an analytic 
circle map. In a higher dimensional dissipative system, xq can be searched either 
as the point visited most rarely by the quasiperiodic orbit or as the limit of points 
Xn, n ^ oo, such that h'^{x) has a minimum at Xn- It could as well be stated that 
h'nixn) ~ Jn{xn)- For a dissipative system, J^(x) ^ as n ^ oo so that a zero 
slope is indeed observed but only considering the limit of an infinitely high iterate 
of the original circle map. 

Furthermore, the calculation shows that the factor of proportion between /i^(xo) 
and Jn{xo) tends to a universal constant as the Jacobian determinant approaches 
zero. In this limit, 

— -- — - ^ a = , , ,„ ,„, , as n — !► oo (1) 

J„(xo) [r7(0)]V'(0) ^ ' 

where ri{x) is one of the components of the universal fixed point pair (^, r]) for the 



standard renormalisation operator T(^, ry) = a{'r], r]o^)a ^ for analytic circle maps 
[8]. MacKay's [12] expansion for r]{x) leads to the numerical estimate a ~ 0.435625. 

The starting point is Bohr's [7] formula relating the derivatives of the first and 
the second iterate of the reduced circle map. Consider a two-dimensional map 
G{x,y) = {gi{x,y),g2{x,y)) with an invariant circle y = c{x) {G{x + l,y) = 
G{x,y) + (1,0)). The original map can be related to a one-dimensional circle 
map by h{x) = gi{x,c{x)) and c{h{x)) = g2{x,c{x)). Differentiating these two 
equations with respect to x yields, after some manipulation, 

^^'^^^ - [gniHx)) + W:^Mi2(^j^,(^^ _ 9i2{h(x)) j^^^ ^2) 



where 



dx gi2{x) gi2{x) 



dgi{x,y) dgi{x,y) 

9ii{x) = — ^ , gi2{x) = — ^ 

dx \y=c{x) dy \y=c{x) 



J{x) = gii{x)g22{x) - gi2{x)g2i{x) 
It is important that the equation for the invariant circle does not appear in (2). 

Eq. (2) can be written in the form 

h'.^{x) = C2{x)h'2{x) + d2{x) (3) 

I introduce a decimation technique by which one can generate from (3) and the 
trivial equation h'2{x) = lh'i{x) + a sequence of equations 

/i^+i(x) = Cn{x)h'^{x) + dn{x) (4) 

with increasing n. Assume that (4) is known for n = i and n = z — 1. Write 

h[j^i{hi{x))h'i{x) = Ci{hi{x))hi{hi{x))h[{x) + di{hi{x))h[{x) (5) 



and split h[(hi{x)) further: 

h[{h,{x)) = c,-i{h,{x))h[_^{h,{x)) + d,-i{h,{x)) (6) 

Replacing h'^_^{hi{x))h[{x) by h'^j^^{x) and using the fact that 

,,, , K^^{x) -d,{x) 

hi{x) = --^ 

c,{x) 

leads finally to 

K+2{x) = Ci+i{x)h'^^^{x) +di+i{x) 

di+i{x) 



with 



Ci+i(a;) = Ci{hi{x))ci-i{hi{x)) 
di+i{x) = 



di{x) 
di{x)[di{hi{x)) + Ci{hi{x))di-i{hi{x))] 



(7) 



Ci{x) 



These recursion relations help in determining the leading asymptotic behaviour 
of Cn{x) and dn{x) as n tends to infinity. First, it can be inductively argued that 
dn{x) ~ Jn{x) (except for n = 1). Recall that di{x) = and d2{x) ~ J{x) 
so that (7) implies ds{x) ~ J{x)J{h{x)) = Js{x). At each level n > 3 of the 
recursion, the leading term in the Jacobian determinant arises from the product 
~ dn{x)dn-i{hn{x)) ~ Jn{x)Jn-i{hn{x))= Jn+i{x). One cau uow proceed to 
determine the asymptotic behaviour of Cn{x) as n ^ oo. Eqs. (5-7) give 



/i^+i(/in(2^)) = Cn+i{x)h'^_^{h^{x)) + ""/] ^^ [h'^_^{hn{x)) - c^{x)] 



dn{x) 



Because dn+i{x)/dn{x) -^ 0, all one needs to know is the asymptotic behaviour of 
h'^_^_i{hn{x)) and h'^_i{hn{x)). I consider here only the critical case with x = xq. 
On the basis of the renormalisation theory [8,9,13], it is expected that 

hn{xo + z)- F„_i - xo ~ a~''r]{a'^z) (8) 
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where the approximation improves with increasing n. This imphes h'^_^_i{hn{xo)) 
-^ r)'{^{0)) = a^ and /i^_]^(/in(a^o)) -^ ^'(^(0)) = '^ where the derivatives have 
been calculated from the fixed point equation. Thus, c^(xo) ^ a^ ~ 1.66 as n 
tends to infinity. 

The asymptotic behaviour of /i^(xo) is solely determined by those of c„(xo) and 
(i„(xo). If h'^{xQ) approached zero slower than (i^(xo), there would exist an A^ such 
that for all n > N , h'^j^i{xQ) > Ch'^{xQ) with C > 1. In this case, h'^{xo) would 
actually keep growing without any limit as n ^ oo which would be contradictory 
to the tendency of the renormalised circle map to approach the universal function 
rj{x). On the other hand, (4) implies that h'^{xo) cannot decay to zero faster than 
dn{xo). In other words, 

h'ni^o) = V^ + 0{Jn+l{xo)) ~ Jn{xo) 

CniXo) 

It turns out to be possible to work out the limit of the factor 

/ \ ^n\X J 

Cn \X jJn\X j 

at x = xo approaching the case in which the Jacobian determinant vanishes. Note 
first that 



gi2ix)c2{x) gi2ix)c2{x)c3{x) 

Eq. (7) implies a recursion relation for enix), n = 3,4,... which becomes very 
simple in the zero-Jacobian-determinant limit: 

en+i{x) ^ en{x)en-i{hn{x)) 

Because C3{x) can be replaced by C2{h2{x))ci{h2{x)) = C2{h2{x)), it is easy to 
write down the form of a general en{x): 

9l2{hn{x)) 1 



en{x) 



912{X) nf=o"' C2{h^(^)) nf=o"' C2(/i'(x)) 



Here I have used the fact that hn{x) mod 1^ x as n ^ oo for the inverse golden 
mean rotation number. Leaving the d-term proportional to the Jacobian out of 
(4), one obtains 

n C2{h\x)) = h'^{h{x)) 

i=0 

where the derivative can be calculated using (8): 

, h'{xo + a-'^r]{a''z))r]'{a'^z) 

^-^'^^"^ + "^^ " JTi^^^T^) 

r]{z) has a cubic critical point at ^ = [8]. Furthermore, also h{xo + z) developes 
such a point in the zero-Jacobian-determinant limit. Expanding all the derivatives 
around z = and letting z ^ leads to 

h^{h{xo)) -^ 

We have thus derived Eq. (1). 

Table 1 shows h'^{xo)/Jn for the dissipative standard map 

k 
gi{x , y) = X + Q + by — —— sin(27rx) 

ZTT 

k 
g2{x, y) = n + by- —- sin(27rx) 

with the constant Jacobian determinant b = 0.5. The critical parameter values for 
the breakup of the "golden" invariant circle can be determined by a dissipative 
version of Greene's residue criterion [14]. xq, j/o is taken as the point where the 
approximating periodic orbits have the largest gap. The calculation of h'^{xo)/Jn 
using the forward recursion relation (4) would be extremely sensitive to the choice 
of the value of /i2(a^o) = h'{xo). An error e in h' [xq] would give rise to an error 
^nr=2 Ci(iCo) in ^n(^o) which would be of the order e 1.66"~^. As J^ decays to 
zero very fast with increasing n (J15 ~ 10~^^^), the error in the ratio h'^{xQ)/Jn 
would soon become astronomical. A better way to calculate this ratio is to apply 



(4) backwards beginning with the approximation h'^{xQ) ~ for some large A^. 
The initial error is very small (~ Jn) and the error in h'^{xQ) is reduced by a factor 
around 1.66 at each step. In fact, a very good estimate for the ratio h'^{xo)/Jn is 
obtained already for n = A^ — 1 if A^ is not very small. The calculation of Cn{xo) 
and dn{xo) by (7) appears to be numerically very stable. Thus, all the error arises 
essentially from the inaccuracy in determining the critical parameter values and 
the point xq, yo- 

Table 1 shows no deviation from (1) although the system is quite far from the 
zero-Jacobian-determinant limit. This could be taken as a hint that Eq. (1) 
would be valid more generically than the derivation would reveal. It would be 
intriguing to see this tested experimentally. If the experimental data enabled one 
to construct the one-dimensional return map, it would be quite easy to calculate 
derivates of higher iterates of this return map by using finite differences. Usually 
the point xq and the Jacobian determinant would not be known. It would be best 
to estimate the smallest slope of each Fibonacci iterate of the reduced circle map 
and calculate h'^h'^_i/h'^j^i. If the Jacobian determinant varied only little along 
the invariant circle, this ratio could be close to the universal constant a. It is 
clear that noise would prevent one from carrying out the calculation for high n. 
Nevertheless, Table 1 suggests that even the lowest order estimate could give a 
reasonable result. 

As to other rotation numbers, I would expect the effective Jacobian again to play 
an important role [15]. The transition to chaos should be observed by monitoring 
the smallest slope of the higher iterate /i*^" (x) of the reduced circle map, where 
Qn would be the denominator of the nth truncation of the continued fraction 
expansion for the rotation number. In the critical case, one would expect this 
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slope to tend to zero as n ^ oo whereas in the subcritical region the asymptotic 
slope should be unity [8,9]. 

Eqs. (4) and (7) can be used to study the conservative case J{x) = 1 as 
well. Both Cn{xo) and dn{xo) have universal non- vanishing limits {xq now 
corresponds to a point on a dominant symmetry line [12]): Coo{xo) ~ 2.1676633 
and doo{xo) ~ —0.4916138. Coo{xo) appears to be determined by the ratio of the 
universal phase space scaling constants [12]. h'^{xo) has a universal positive limit 
doo{xo)/{l — Coo(a^o)) ~ 0.4210236. The estimate obtained by setting J^ = 1 in 
(1) deviates only about 3 % from this true value. 

The decimation technique introduced in this letter is readily applicable also to 
other problems, e.g. to the discrete quasiperiodic Schrodinger equation [16]. Eq. 
(3) can be interpreted as a discrete eigenvalue equation with h'^ representing 
the wave vector i/jp^ at site Fn [17]. The modulating potential is included in 
C2- Assuming the normalisation condition •i/'o = 1 one can take d2 = —1. The 
present approach is appropriate when the potential has the frequency ( relative to 
the underlying lattice. Infinite products of transfer matrices usually diverge [18] 
whereas, by writing down recursion relations similar to (7), it is possible to find a 
bounded limiting behaviour for the coefficients Cn and dn [17]. 

The fact that the slope of the reduced circle map depends on the effective Jacobian 
is in a nice agreement with the conjectured mechanism for the breakup of having 
a tangency between the invariant circle and its stable foliation [8]. It is natural 
to think the contraction on the stable foliation to be proportional to the effective 
Jacobian determinant. At the point of tangency, the slope of the associated circle 
map should therefore be proportional to the Jacobian. 
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Table 1. Subsequent estimates for the factor h'^{xo)/Jn resulting from the 
critical dynamics of the dissipative standard map (6 = 0.5, k — 0.9788377790, 
O = 0.3058769514). 



n 


h'r,{xo)/Jn 


2 


0.5280 


3 


0.4333 


4 


0.4548 


5 


0.4146 


6 


0.4384 


7 


0.4191 


8 


0.4385 


9 


0.4272 


10 


0.4387 


11 


0.4312 


12 


0.4380 


13 


0.4331 


14 


0.4373 


15 


0.4339 
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